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already been successfully used in a number of astrophysical research projects. 

Index Terms — 1.3. 5.f Modeling packages, 1.4.8 Scene Analysis, J.2.c Astronomy, J. 2.1 Physics 

> 



1 Introduction 

The interpretation of astrophysical data often depends 
strongly on the knowledge of depth information along 
the line of sight. In most cases, however, this is the least 
well known information. That is true for the distance 
and especially for the position of substructure within 
an object. The development of effective methods for 
the reconstruction of the 3D structure of astrophysical 
objects is therefore an issue of growing importance in 
astronomy. Photographic images only provide a two- 
dimensional integration of the emission and absorption 
along the line of sight. The depth information is therefore 
flattened. Sometimes, symmetry properties combined 
with a favorable orientation of an object provide suffi- 
cient information to visually deduce what the structure 
must be. This can be the case for planetary nebulae 
and has been used to automatically reconstruct the 3D 
structure (Leahy ||l4l|, Magnor et al. |^, |^, Lin^u 
et al. 1 16 1, |17|). If no such symmetries are present, then 
the depth information must come from other types of 
information, which usually depend on a fundamental 
physical model for the object class that is considered. 
This information could be the velocity field, e.g. in a 
radially expanding nebula a mapping between velocity 
and position exists. However, for some objects - such 
as turbulent interstellar clouds - such a mapping is not 
possible. 

Much of observational astrophysics research involves 
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physical modeling with limited constraints to deduce 
physical properties of the observed objects. Astrophysi- 
cists measure a limited number of physical properties of 
an object, via electromagnetic waves, to which a physical 
model of the phenomenon is then fitted. Such models are 
usually not unique. 

Most of the astrophysical modeling effort tends to 
gravitate towards massive parallel supercomputing for 
dynamical simulations. Analysis and visualization of 
such simulations are done separately and often are 
complex and computationally intense processes them- 
selves ITT) . While such simulations produce insight into 
generic astrophysical processes, they are rarely suitable 
for elucidating the properties and structure of particular 
objects. 

Knowing the properties of individual objects is essen- 
tial when a single object class, e.g. planetary nebulae, 
shows a large variety of presentations in images. A 
serious difficulty for the modeling of particular objects 
is our fixed vantage point on Earth which restricts all 
observations to be along a single direction (up to the 
parallax provided by the Earth's orbit, which is negligi- 
ble for typical distances to astronomical nebulae). This 
is in strong contrast to, e.g., medical imaging where 3D 
information is recovered from observations from multi- 
ple directions around the subject. For any astronomical 
object beyond the solar system, we are able to observe 
only one 2D projection of its actual 3D volumetric shape. 
For the correct physical interpretation of observational 
data, information about the object's actual 3D shape has 
to be available ||3|, (15), |18|, |24|, |34|. Obtaining new 
structural information and insight on particular objects 
is the main purpose of the application that we present 
in this paper. Our approach to modeling individual 
objects is very different from previous methods. The 
application that we present (called Shape) becomes essen- 
tial when automatic reconstruction methods fail because 
theoretical or observational constraints are insufficient. A 
lack of constraints for an automatic reconstruction is at 
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least partially compensated by scientific user judgement. 
Often the available constraints are sufficient to test one 
or more hypotheses about the structure of an object. 

In such cases, rather than a reconstructive, a construc- 
tive morpho-kinematical modeling approach is more 
suitable to reconstruct the structure and velocity field. 
The term morpho-kinematical is applied to modeling 
that involves only structural (morphological) and veloc- 
ity (kinematic) information. This is in contrast to dy- 
namical simulations, which include the effects of forces 
and temporal evolution from a set of simpler initial 
and boundary conditions. In general, the outcome of 
dynamical simulations is not predictable in detail and 
very hard to tune to a specific object. 

Conventional morpho-kinematical modeling uses 
hard-coded mathematical descriptions of the objects, 
processes and boundary conditions. Therefore, the user 
needs at least basic programming skills in the particular 
language of the code. Modern 3D modeling software of 
the graphics industry shows that such modeling can be 
done effectively without user programming intervention. 
Although such software can visualize gas-dynamical 
processes, it is inefficient and its usefulness for astro- 
physical processes is limited. The general workflow of 
such systems seemed, however, very suitable for mod- 
eling particular astrophysical objects [33]. 

Following the technique of modern interactive 3D 
graphics systems, we have developed Shape with special- 
ized functionality for interactive astrophysical modeling 
on single desktop or laptop computers. The primary 
purpose of Shape is to interactively generate 3D models. 
However, in contrast to conventional astrophysical mod- 
eling tools, it integrates the visualization and analysis 
of the model into the same system. Direct access to the 
model data at any stage of the modeling process allows 
for effective comparison with the observed data in a 
variety of ways within the feedback loop of the iterative 
workflow (Figure 

In artistic work on astronomical topics, commercial 
tools like Maya or 3D Studio Max are frequently being 
employed |4|, |22|, |37|. Professional animation tools 
are designed to assist in creating realistic 3D scenes 
of familiar environments. Unfortunately, when used for 
scientific work they display serious shortcomings. Espe- 
cially volume rendering with mesh structures and parti- 
cle systems are very different from the physical correct 
radiation transfer needed for reliable interpretation of 
astrophysical phenomena. A qualitative and quantitative 
comparison of such models with real objects is not possi- 
ble. The key problem preventing their use in astrophysics 
research is the inability to produce the type of render- 
ings that are comparable to the observations obtained 
with telescopes and other scientific instrumentation like 
spectrographs (e.g., for Doppler-shift measurements). 

With Shape we remedy most of the shortcomings of 
previous astrophysical reconstruction systems by apply- 
ing the powerful structure modeling techniques of com- 
mercial animation suites, while adding the information 
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Fig. 1: The Shape modeling tool that we present here 
provides a unified framework for the whole modeling 
and visualization workflow. The user employs his phys- 
ical knowledge to construct an initial model which can 
be visualized and compared to observational data in 
several ways, allowing for easy interactive and iterative 
refinement of the model. When all necessary physical 
information is reflected in the model, its parameters can 
be automatically optimized, minimizing the difference 
between the model and the observational data. The final 
model can then be used to generate various types of 
graphical output. 



output and processing systems that are necessary for 
astrophysical research applications. We go beyond the 
current commercial rendering techniques by using phys- 
ically more accurate modeling of the radiation transfer 
from the sources to the observer. 

In this paper we first comment on previous related 
work in Section |2j on the type of observational data that 
are used for this work in Section |3] and then introduce 
the Shape system. In Section |5] we show three examples 
of previously published research applications of Shape, 
before giving an outlook on future developments and 
our conclusions in Sections |6] and |7j respectively. 

2 Related Work 

Many different approaches have historically been used 
to simplify the modeling of astrophysical objects. Be- 
sides the general reconstruction techniques applying 
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to a wider range of volumetric objects (cf. flS]), spe- 
cialized approaches have been developed that exploit 
the peculiarities of the astrophysical case. Among them 
are automatic methods that make use of Doppler shift 
measurements (cf. the section on observational data) or 
symmetry assumptions as well as user-driven modeling 
systems. 

2.1 Automatic Reconstruction Methods 

Many extended astrophysical objects show a strong cor- 
relation between the velocity and position of the emitting 
gas relative to some local reference point, which often 
is the center of a star or stellar remnant. For example, 
material that has been ejected ballistically with different 
velocities from the same source (e.g. by an explosion) 
will naturally evolve such that - after a time period that 
is long enough compared to the duration of the ejection 
process - the faster regions have moved farther away 
from the source and an approximately linear dependence 
between velocity and distance from the source is estab- 
lished (see Figure |2j top). 

Doppler shift methods like the one by Sabbadin et 
al. 1 27 1, 1 28 1 make use of this correlation in order to 
derive depth information from Doppler-shift data. If the 
assumption of linear dependence between the position 
and velocity vector holds, a linear mapping exists be- 
tween the Doppler-shift (i.e. velocity along the line of 
sight) and the position along the same direction (see 
Figure [2j bottom). In this case the resulting models 
are accurate within the limits of the accuracy of the 
Earth-bound observational data. Unfortunately, many 
objects contain several different kinematic subsystems 
which may have different relations between velocity and 
position. Some also show complex interactions with their 
local environment which may further complicate the 
velocity law |32|. Furthermore, these methods require 
an almost complete coverage of the object with regularly 
spaced observations of the Doppler-shift, which require 
special observing programs. Such homogeneous data 
sets are rarely available. 

Other algorithms that are more oriented towards high 
visual quality than physical accuracy of the results are 
based on symmetry constraints (Magnor et al. |19|, |20|, 
Lintu et al. |16|, |17|, Wenger et al. |35|). Many astro- 
nomical nebulae show an inherent spherical or axial 
symmetry due to their evolution from more or less 
symmetrical sources. This symmetry assumption may be 
used to reconstruct the missing spatial dimension |14|. 

In many cases, however, these simple symmetries are 
disturbed by statistical effects or external influences, and 
more complex symmetries as well as turbulent structures 
arise which cannot be modeled in such a generic way. 

2.2 User-Driven Systems 

The above-mentioned shortcomings of fully automatic 
reconstruction approaches may be avoided by resorting 




Fig. 2: The top panel shows the outline of a bipolar neb- 
ula with position and velocity vectors marked. Most 3D 
reconstruction systems for expanding astrophysical neb- 
ulae are based on the assumption that the velocity vector 
is proportional to the position vector. The Doppler-shift 
of a spectral line due to the velocity component along the 
line of sight (vz, bottom panel) provides a direct mapping 
to the position along the line of sight (z). The unknown 
constant of proportionality can often be determined from 
additional information, e.g. symmetries that are inherent 
in the object structure. 



to interactive modeling techniques. The interactive cre- 
ation of a model aiming to reproduce a given single 
image is a common task, but most existing solutions 
are not well suited for the modeling of emissive trans- 
parent objects that are prevalent in astronomy and do 
not allow for representation of velocity information and 
spectral data. Among the tools that most closely reflect 
our modeling approach are the interactive approaches 
of Debevec et al. |5|, Francois and Medioni |8| and 
Zhang et al. |38|, all of which expect some kind of user- 
specified coarse geometry or a set of user-defined geom- 
etry constraints which is then automatically converted 
into a full three-dimensional model that best fits the 
provided image under the given constraints. The idea of 
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Fig. 3: Panel "a'' shows the projected synthetic image 
of a bipolar nebula comparable to the one outlined in 
Figure ^ with a similar orientation to the observer. The 
color coding is according to the velocity along the line 
of sight (blue is approaching and red is receding). The 
marked synthetic spectrograph slit was used to generate 
the P-V diagram in panel 'V\ Since the velocity field 
has been assumed homologous, the P-V diagram can be 
stretched along the horizontal velocity axis such that its 
outline corresponds to a cut through the object along the 
line of sight (panel "c'O- 



automatic optimization of a parameterized (deformable 
or "morphable'O model has been successfully employed 
in the works of Montagnat and Delingette |23J and 
Romdhani and Vetter |26|. 

Some entirely manual modeling aids have been de- 
veloped for the astrophysical use case, in which an 
astronomer defines a model representing his or her 
high-level knowledge about the object in question. This 
process will usually involve iterations of modeling, com- 
parison of the rendered model to actual observational 
data (especially spectral data), and refinement of the 
model, until all observational facts can be explained by 
the physically plausible model. 

An early tool for the rendering part of the astrophys- 
ical modeling process was our earlier work [llj which 
was able to reproduce many standard forms of observa- 
tional data from a given model. Similar codes have been 
used by Santander-Garcia et al. |31| and Hajian et al. |10|. 
The model itself, however, still had to be hard-coded 
into the program, making the modeling part inherently 
cumbersome. Steffen and Lopez |33| later incorporated 
their spectral renderer into a commercial modeling sys- 
tem as a plugin. This simplified the modeling process 
to a large extent, but performance and usability were 
still far from the quality of an integrated modeling and 
rendering system. 

The only astrophysics-suited tool that we know of 
which follows a paradigm similar to that of our earlier 
work is Hydra which is under development at MIT ||6|, 
ijT^. Contrary to Shape, Hydra focuses on X-ray data 
instead of the near-visible wavelengths which are suit- 
able for realistic visualizations aimed at public media 
like planetaria |21|. Also, the tool does not provide an 
interactive modeling system, but models are defined 
using scripting and Constructive Solid Geometry (CSC). 



3 Observational Data in Astronomy 

Many astrophysical objects (e.g. stars or emission nebu- 
lae) contain large amounts of ionized gas. This ionized 
gas (or plasma) emits light at a number of well-known 
wavelengths, depending on the chemical composition of 
the plasma. Using filters for selected wavelengths, two- 
dimensional projections of the spatial distributions of 
different ions within such an object can be measured 
as photographs. While the spatial resolution of earth- 
bound measurement devices is limited by the perturbing 
effects of the atmosphere, the advent of space telescopes 
such as the Hubble Space Telescope (HST) made high 
resolution imagery with high sensitivity available for a 
large number of objects. 

Another key source of astronomical data for 3D recon- 
structions are Doppler-shift measurements. The underlying 
physical phenomenon is the Doppler-effect which causes 
the wavelength of light to change when the emitter and 
the receiver have a nonzero relative velocity towards 
each other. With high-resolution spectrographs, it is pos- 
sible to measure the shift of any observed spectral line 
with respect to its known reference wavelength, and 
therefore the relative velocity of the emitter along the line 
of sight. Because many complex astrophysical objects 
also have complex internal velocity fields, the resulting 
distortions of the spectral lines are an important addi- 
tional source of information. 

To capture spatially varying wavelength information, 
usually a narrow slit is used to select only a narrow, but 
long region of the object. The light is then dispersed in 
wavelength along the direction perpendicular to the slit. 
The result is recorded on an image detector and leads 
to a position-velocity (P-V) diagram. The image intensity 
is the object intensity as a function of position along the 
slit and wavelength (Figure |3|. 

4 The Shape System 

The motivation for the development of Shape was to 
be able to reconstruct the 3D structure of astrophysical 
objects based on data sets that do not allow an automatic 
reconstruction. This may be because of a lack of suffi- 
ciently detailed data or because the basic assumptions 
for the reconstruction algorithm are not met by the 
object. Therefore, a user-driven interactive approach was 
adopted. The minimum general requirements that such 
a system should have are those present in previous 
systems, which basically are: 

• Tools to define a spatial emissivity and velocity field. 

• Establish a mapping between Doppler-shift and 
position. 

• Produce output that can be compared with observed 
images and spectra. 

• Shape: The tools to define the spatial structure and 
velocity field should be interactive 

The first three requirements are for previously existing 
systems. The last item, in italics, is the fundamental 
additional distinctive requirement for Shape. 
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Fig. 4: The Shape graphical user interface contains several modules dedicated to different tasks. Here, the interactive 
3D modeling window, the screen for rendering and comparison to observational data and the control panel for 
automatic optimization are shown, respectively, from left to right. 



To fulfill these four requirements, three major Shape 
modules are typically used when constructing a model to 
fit given observational data. The workflow is illustrated 
in Figure First, the interactive modeling environment is 
used to specify the user's assumptions about the object's 
geometry, velocity and emissivity. This model is then 
input to the rendering module which emulates how the 
modeled object would be observed through a physical 
measurement device such as a telescope with a spec- 
trograph. Visualizations have been designed such that 
they are an aid to obtain new physical insight and /or 
be comparable to actual astronomical data. The user then 
interactively refines the model until it qualitatively fits 
the observational data. In a final automatic optimization 
step, the numerical parameters of the model are varied 
to also quantitatively fit the model in the best possible 
way. 

The interactive modeling approach makes it important 
to keep the system highly functional on single comput- 
ers. Parallelization therefore focuses on the application 
of multi-core threading and, in the future, on parallel 
simulation and rendering on graphics processing units 
(CPUs). 

4.1 Java Implementation 

For the choice of Java as the programming language for 
Shape the following three initial criteria were decisive. 
Astrophysicists use a range of operating systems, mainly 
Linux, Mac OS X and Windows. Therefore, first, portabil- 
ity of the software and file formats was essential. Second, 
since the software was to be developed over a consid- 
erable time and updated regularly, Java Webstart seemed 
an excellent way to handle user-friendly updates. Third, 
all essential Application Programming Interfaces (APIs) 
are readily available. 

In practice, however, there are a number of compatibil- 
ity issues, especially with Mac OS X. The necessary Java 
software has become available only recently and only for 
64-bit systems. Minor non-essential compatibility issues 



remain with Linux and 64-bit Mac OS X systems, which 
can be expected to disappear as these operating systems 
and Shape evolve. 

Shape requires four basic Java components: the Java 
Virtual Machine (JVM), Java3D, Java Advanced Imaging 
{J AT) and Java Media Framework (JMF). The 3D modeling 
environment is implemented with JavaSD, and JAl has 
been used for some image processing tasks. JMF allows 
the encoding of animated image sequences into common 
movie formats. 

Depending on the selected algorithm and the com- 
plexity of a model, simulation times may range between 
fractions of a second and about one or two minutes on 
a current Intel Core Duo processor with 2 GB RAM, 
but they usually stay below a few seconds for the most 
common cases with typical 256 pixel image sizes. 

4.2 Interactive Modeling 

The modeling system that we present here has to solve 
two key problems that are common in existing astro- 
nomical 3D reconstruction tools. Since the spatial and 
spectral resolution of current observational techniques 
reveal a great wealth of irregular detail in gaseous neb- 
ulae, simple mathematical descriptions of the structures 
have become impractical as the number of parameters to 
handle can go into the hundreds. We want our system 
to be able to model the complex multipolar, knotty or 
filamentary structures that are commonly observed in 
gaseous nebulae. Also, most earlier solutions have been 
based on "off-line'' modeling of the object being recon- 
structed, followed by an entirely separated rendering 
step. We want to allow for an interactive model-render- 
compare feedback loop in order to speed up and simplify 
the modeling process. 

These problems lead us to the specification of a mesh- 
based interactive construction software similar to com- 
mercial modeling tools. The main difference between 
our tool and conventional 3D animation packages is 
that Shape produces not only images but also spectral 
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information which is presented in a way to be directly 
comparable to various types of astronomical observa- 
tions, and that its modeling environment focuses on the 
structures and symmetries that are commonly present in 
astronomical objects. 

The Shape program consists of two main components: 
A 3D modeling view in which geometry and behavior 
of the model are defined, and a 2D screen where the 
simulated appearance of the model can be compared to 
observational data (Figure |4|. 

The model properties are defined basically in three 
steps. First a structural selection of regions in the spatial 
domain is set using mesh surfaces and volumes. Multiple 
meshes can be used to model complex structures. Physi- 
cal properties are assigned as a function of position in the 
spatial domain. Different functions can be used for dif- 
ferent sub-structures, i.e. meshes. A sampling process in 
the selected spatial mesh regions then reads the physical 
properties. The sampling can either be done at particle 
locations on the surface or within the volume of the mesh 
or at the locations of the voxels of a regular 3D grid. 
Sampling particles can be distributed randomly on the 
surface or the volume of a mesh. The particles can also 
be used to visualize the vector field, attaching to them 
a velocity vector. Although these particles have physical 
attributes, they should not be confused with physical 
atomic particles. Here they serve only the purpose of 
visualization and of defining a position for the sampling 
process. 

For the mesh creation process, a number of suitable 
primitives (sphere, torus, cone, cube, etc.) are available 
which are then deformed using a selection of modifiers; 
importing meshes from other software is also possible. 
The primitives can be used either as volumetric objects, 
as infinitely thin shells, or as volumetric shells with user- 
specified thickness. In addition to scaling, translation 
and rotation, the available modifiers include operations 
that are parameterized along a given axis (the squeeze, 
squish, shear and twist modifiers). Boolean combina- 
tion of different primitives allows for constructive solid 
geometry modeling. Additionally, mesh vertices can be 
adjusted manually either individually or in groups. 

Among others, the spatial distribution of density, color 
and velocity can be conveniently defined as a function of 
position. The user can choose among a set of predefined 
common spatial dependencies with adjustable parame- 
ters or create custom functions in different coordinate 
systems (cartesian, cylindrical or spherical). These func- 
tions are assumed to be separable into functions of the 
chosen coordinates, e.g. /(r, 6>, ^) = fr{r) • fe{0) • /^{(j)). 

The different coordinate functions (/r, fe and Z^, 
in that example) can be defined either analytically or 
as piecewise linear functions that may be graphically 
edited (like animation curves in conventional animation 
software packages). For the analytic functions, the user 
can define the function by typing in a formula. The 
formula is fed into a mathematical text interpreter and 
is calculated and displayed immediately as a graph. 



Predefined functions for, e.g., the velocity field may 
include common movements such as radial expansion, 
gaseous disk rotation and solid body rotation, and ran- 
dom distributions. More complex functions can be con- 
structed by sequentially combining them with addition 
or multiplication. 

All modeling and visualization can be done in either 
arbitrary units or in various actual physical units that 
are suited to the problem, e.g. the apparent size of 
the object in arcseconds, sizes and distances in parsec 
or astronomical units, velocities in km/s and the like. 
The modeling interface not only displays the mesh and 
particles, but also the particle velocities in the form 
of vectors color-coding their velocity along the line of 
sight, which is responsible for the observed Doppler- 
shift. The projection of these vectors on the sky is directly 
observable in some objects as tangential motion, after 
comparing observations with a sufficiently large time 
interval in between. This provides further kinematical 
constraints (see also Section [5^ Figure |9|. 

Since the particle velocities can be specified in a model, 
the time evolution of the object may be predicted, assum- 
ing ballistic expansion (i.e. constant velocities). A time 
modifier calculates the future or past positions of the 
particles for a given time interval. Using this feature, the 
age of a nebula and its short-term structural evolution 
can be estimated. 

Finally, our tool also allows importing data from 
external simulations in order to visualize and analyze 
the results (cf. the examples section). This makes it also 
possible to use other software for specialized modeling 
tasks (e.g. for the generation of sophisticated noise dis- 
tributions). Also, the models created with Shape may be 
exported and used as input for external simulation or 
visualization software. 

4.3 Image Rendering 

Several renderers are available, which serve different 
purposes according to the type of object that is modeled 
and the adopted workflow. Many astrophysical objects 
are optically thin, i.e. transparent. This fact has been 
used in three renderers {particle, grid & mesh renderers). 
They provide faster rendering than the physical renderer, 
which takes into account opacity and other radiation 
transport effects. The physical renderer is still experimen- 
tal and under development. Its details will be described 
elsewhere. 

The particle renderer uses a random particle distribu- 
tion to sample the model emissivity and velocity space. 
The values are added directly to the image pixels or 
P-V diagrams according to their projected position and 
velocity along the line of sight. 

In the grid and mesh renderers, a regular grid is set up 
in world space that is aligned with the line of sight. In 
the grid renderer, particle positions are used to sample 
the physical properties of the object. The density of 
the particles are distributed in the 3D grid. The voxels 
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Fig. 5: The image on the left is a rendering of a 
hydrodynamical simulation after being imported into 
Shape. The middle shows the Shape mesh fitted to this 
image and a small fraction of the sampling particles. The 
panel on the right is the rendering of the Shape model of 
the hydrodynamical simulation after manually fitting an 
emissivity distribution in cylindrical coordinates ||32J. 

are then rendered and/or output for external purposes. 
For the mesh renderer, the sampling of the physical 
properties is done by searching for the voxels of the 
grid that overlap with the model's mesh. If the mesh 
is a volume, then the center point of the voxel is used as 
the sample point and emission is calculated from the 
fraction of the voxel that is located within the mesh. 
To reduce aliasing effects, the position can be jittered 
within the limits of the voxel. If the mesh is a surface, 
the fraction of each mesh segment in a particular voxel is 
determined. Assuming a small but finite thickness of the 
surface, the fraction of the volume occupied in the voxel, 
and hence the emission, can be determined. The emission 
from each sub-object is added to the grid and rendered 
separately (no mixing is done). Finally, the emission from 
each voxel is projected onto the image plane. 

In the images, the emission is integrated along the line 
of sight (z-axis) regardless of their velocity, and the re- 
sults are plotted in the image plane (xy-plane). Position- 
velocity diagrams (P-V diagrams), however, only take 
emission in a given x range (within the spectrograph 
slit) into account. The intrinsic spectral line width is 
assumed small compared to the resolution of the P- 
V diagram. The y coordinate is the position along the 
slit. The object's emission is distributed according to its 
position along the slit and the velocity component along 
the line of sight (v^, see Figures |2] and |3|. 

4.4 Synthetic observations 

An important task in the process of reconstructing an 
object based on observational data is, of course, the com- 
parison between these data and the model. In order to 
make this comparison reliable, the rendering algorithm 




Fig. 6: Position-velocity diagrams are shown for the 
hydrodynamical simulation (grey scale) and Shape model 
(contours) of Figure (|5| for three different viewing an- 
gles. Note that the outline changes with the viewing 
angle. Furthermore, the outline of the first P-V diagram, 
which corresponds to the viewing angle used in the 
rendering of Figure (|5|, is different from that of the 
image. The middle P-V diagram shows point-symmetry, 
whereas the others are mirror-symmetric. Both observa- 
tions are indications for deviations from a homologous 
expansion of the nebula |32|. This shows that reconstruc- 
tion methods that rely on homologous expansion will 
yield qualitatively inaccurate results for such a nebula. 



has to take into account not only the model itself, but 
also the properties of the measurement devices. 

For comparison with observations. Shape includes 
three main 2D data types that are rendered from the 3D 
model: images, P-V diagrams and channel maps. For 
the images and P-V diagrams, several coloring schemes 
are available that display different types of secondary 
information. First, the color of a substructure can be used 
as assigned in the 3D model. This helps to distinguish 
the different sub-systems of the object, especially in the 
spectra. Second, the velocity along the line of sight as 
given by the Doppler-shift can be color-coded as simple 
red /blue to distinguish regions that move towards or 
away from the observer, respectively. Another option 
maps the Doppler-shift to a rainbow color range, such 
that a more detailed appreciation of the line-of-sight 
motion can be achieved. 

Figure (|3| shows the rainbow display option for image 
(a) and P-V diagrams (b & c). In panel (a) the image of 
a synthetic bipolar structure is shown as seen with a 
structure and inclination angle similar to the outline in 
the schematic diagram of Figure (|2|. The color-coding 
maps the Doppler-shift to the rainbow colors from red 
to violet, with red receding and blue approaching the 
observer. The width of the spectrograph slit has been 
marked. Panel (b) of the same figure is the resulting P- 
V diagram with a scaling in velocity (horizontal) that 
is typical for high-resolution spectral observations of 
bipolar planetary nebulae. In this object the velocity 
vector is proportional to the position vector (homologous 
expansion), such that for a suitable velocity scaling, the 
P-V diagram should reproduce the object structure along 
the line of sight. This is confirmed in panel (c) where the 
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velocity scale has been stretched accordingly. If there are 
deviations from a homologous expansion, the structure 
of the P-V diagram can not be stretched to yield an 
accurate representation of the object structure along the 
line of sight (see the examples section below). 

A third image type that is common in astronomy is 
the so-called channel map where only a given velocity 
(or wavelength) range is projected onto the x?/-plane; 
Shape is able to produce this output, optionally pre- 
senting maps for many different velocity ranges next 
to each other. This observing mode is typical for radio 
wavelength observations of spectral lines of molecules in 
dense star forming regions and young planetary nebulae. 

The channel map mode in Shape can also be used 
to display light echoes, where each image is a channel 
of equal light travel time. Such light echoes may occur 
when a star experiences a sudden and short outburst. 
The light of the outburst may then be reflected by 
surrounding dust. The appearance of the nebula changes 
on a time scale corresponding to the time it takes to 
propagate from the star through the surrounding dust. 
At a given time, we on Earth can only observe the radi- 
ation that has taken the same amount of time from the 
star to the nebula and then to Earth. Therefore, the light 
echoes represent channels of equal light travel time. The 
geometrical shape of such a light echo in space is that of 
a paraboloid. It has the property of equal distances from 
the position of the star (the focus) to a given position 
on the paraboloid and then to the observer. The channel 
map mode allows one to calculate and display such light 
echo images. A well known example of such a light echo 
is the one observed around the star V838 Monocerotis 
after its large outburst in 2002 |2|. 

Some instrumentation only registers particular spec- 
tral ranges or has a varying sensitivity at different wave- 
lengths. Therefore a set of filters can be applied during 
the rendering process, that filters the emission according 
to a customizable set of rendering filters. These adjust or 
exclude emission according to the specifications of the 
filter (see below, the example of nova RS Ophiuchi). 

Final post-processing steps on the rendered raw image 
and spectral information are used to model the prop- 
erties of observational devices, e.g. the influence of the 
atmosphere or the aperture. The most important of the 
atmospheric and instrumental effects are the spatial and 
spectral resolution. These are modeled by applying gaus- 
sian convolution to the spatial and spectral dimensions. 
Effects like complex point-spread functions and ghost 
images from reflections inside the instrumentation are 
not taken into account. Other filters like gamma cor- 
rection or logarithmic scaling make high dynamic range 
data manageable. Contour levels of the image intensity 
may also be applied. The rendered image and spectrum 
can then be displayed and interactively blended with the 
observation to visually detect differences. 

When the model satisfies the demands of the user, a 
number of visualizations can be done that focus on the 
presentation of the results instead of being a modeling 



aid. For example, images may be rendered from differ- 
ent viewpoints, and image sequences and videos may 
be generated by changing display parameters (e.g. the 
camera position or the slit position of a P-V diagram). 

4.5 Automatic Optimization 

Shape includes an automatic optimization module which 
minimizes the difference between the rendered model 
and the observational data by tweaking the values of 
observing parameters or those of the modifiers in the 
model. The user can influence this process by selecting 
the parameters to be optimized and by limiting their 
values to a plausible range, as well as by optimizing 
different parameters in order of priority until the model 
adequately fits the observations. 

Since a model may contain hundreds of parameters 
that influence each other, the user is presented with a 
hierarchical list of optimizable parameters from which 
a reasonable subset (usually only a few parameters) 
are then selected for optimization. This interactive pa- 
rameter pre-selection not only speeds up the following 
algorithm, but it may also be seen as an additional 
way for the user to emphasize certain properties of the 
model: Since changing the value of any one of two 
independent parameters might make the model explain 
the observational data equally well, there is ambiguity 
in any optimization process that can only be resolved 
by the user. Specifying ranges of valid parameter values 
is an even more powerful method of user intervention, 
because not selecting a parameter for optimization is 
eventually equivalent to selecting an infinitely small 
value range. Many parameter types presuppose certain 
boundaries for their values. For example, in a squeeze 
modifier which scales the distance of vertices from a 
given axis, the scaling coefficients are necessarily pos- 
itive. For the optimization one can use images, P-V 
diagrams or data plots, of which only one can currently 
be optimized at a time. 

4.6 Automatic reconstruction 

Although Shape has been conceived for interactive recon- 
struction, there is an automatic reconstruction module. 
It is intended to provide a first look at the structure of a 
complex object for which a suitable data set is available. 
From the input data it generates a particle system that 
follows the brightness distribution of the data. 

Our system has several advantages compared to other 
automatic reconstruction systems that are based on a 
mapping between Doppler-shift and position along the 
line of sight. First, it can use two different data types: P- 
V diagrams and channel maps. Second, it is not restricted 
to the assumption of a homologous expansion. It can 
use any radial velocity law as long as it is monotonic 
and provides a one-to-one mapping from Doppler-shift 
to position along the line of sight. 

The input data in the form of P-V diagrams, in con- 
trast to the tomographic method by Sabbadin et al. pO) , 
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are assumed to come from parallel slits. Parallel slit 
settings are much easier to obtain and process, since they 
can cover the complete object without gaps or super- 
positions. Similarly, a set of channel maps can be used 
to recover the 3D structure. Currently, no interpolation 
is performed if there are gaps in the spectrograph slit 
coverage. The gaps can, however, be filled by assum- 
ing that the slits are wider and touch each other. Slit 
positions and widths for the input P-V diagrams can 
be set individually. As usual, the exact scaling between 
Doppler-shift and position must come from some sym- 
metry information about the object or some part of it. It is 
clear that such a reconstruction can only provide a first 
approximation to a detailed interactive reconstruction, 
especially if the assumed velocity law is not accurate for 
any part of the object. 

4.7 Plot, Animation and Movie Modules 

There are three auxiliary modules in Shape that are fully 
integrated in the main modeling process. They are the 
modules for plotting data, animating model parameters 
and the display of animation results. 

The plot module can display graphs of a variety of 
data from the models, including ID spectral line profiles. 
While the spectral line profiles require the rendering of 
the model, other data types can be obtained directly from 
the particles in a 3D model. For instance, when plotting 
velocity along the line of sight as a function of position 
along the spectrograph slit, a real-time preview of the 
P-V diagram is obtained in the form of a scatter plot of 
these particle properties. Such data can be visualized in 
real-time, as the user changes, for instance, the camera 
orientation or the slit position in the other modules. 

The animation module was inspired by the equivalent 
curve editor modules of conventional 3D animation 
software. Here almost all object and camera parameters 
can be animated as a function of time using either 
manual animation curve editing or algebraic functions. 
In addition to educational visualizations, this allows 
the scientist to effectively explore the parameter space 
by producing sequences of outputs by automatically 
changing the parameters in a controlled manner. The an- 
imation module has an interactive time-line that updates 
the object structure in the 3D module in real-time. After 
rendering a sequence, the results can be viewed in the 
movie module. The movie module can load and simul- 
taneously reproduce an arbitrary number of animation 
sequences, which may include images, spectra and plots. 
This is very helpful when analyzing various types of 
outputs as the parameters of a model change. 

5 Results and Example Applications 

In this section we show three examples of Shape models 
that have been published in the astrophysical research 
literature. The first example is an analysis of hydrody- 
namical simulations. These have been used to validate 
Shape and extract information about deviations from a 



homologous expansion that can be expected in planetary 
nebulae. Using Shape these deviations have been reduced 
to a small set of parameters for later use in models of real 
objects, like the one of the planetary nebula NGC 7009, 
which is presented after the hydrodynamical simulation. 
The third example is that of Nova RS Ophiuchi, which 
has been constructed by a user that is independent of 
the Shape developers. It applies features that have been 
implemented upon request by the user. 

For this work. Shape has served as a completely inte- 
grated tool for modeling, simulation, analysis and visu- 
alization, setting it apart from any existing commercial or 
research tool. The results shown in this section represent 
original astrophysical research results which would have 
been very hard or impossible to obtain without Shape. 

5.1 Validation of Shape with Hydrodynamical Simu- 
lations 

In addition to direct modeling of astrophysical obser- 
vations. Shape has been applied to the analysis and 
visualization of morpho-kinematical aspects of numer- 
ical hydrodynamical simulations. Such simulations have 
also served to validate the functionality of Shape, since 
their properties are known in full. Steffen & Garcia- 
Segura |32| have used Shape to characterize the velocity 
field of numerical simulations of some basic types of 
planetary nebulae. 

Parametric descriptions of the velocity field that ex- 
tend commonly assumed velocity fields in morpho- 
kinematical modeling of axisymmetric objects where 
derived. 

For this analysis, the simulated hydrodynamical data 
have been filtered according to density. The densest 
elements correspond to the shell that is usually observed 
as the brightest region in a planetary nebula. They have 
been imported to Shape as a particle system, including 
their velocity information. A 3D mesh and density distri- 
bution was then fitted to the large-scale structure. Figure 
(|5| shows a comparison of the integrated emission mea- 
sure (i.e. density squared) from one of the simulations 
(left) and that of the corresponding Shape model (right). 
The mesh structure and particle distribution (10 % of the 
actually used particles) are shown in the middle panel 
of the same figure. 

The authors allowed the presence of a poloidal ve- 
locity component. This component is perpendicular to 
the radial velocity component and points towards the 
symmetry axis along the "longitude'' of a spherical 
coordinate system. As a function of angle from the 
symmetry axis, the poloidal velocity can be described 
by three linear segments with the condition of zero 
magnitude on the symmetry axis and equator. The radial 
component can also be described by two or three linear 
segments. These descriptions of velocity fields capture 
the kinematical properties of nebulae quite accurately, 
and are simple enough to be easily implemented in 
morpho-kinematical or photoionization codes (Figure l6|. 
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To verify the kinematic modeling functionality of 
Shape, the velocity field of the simulations was manually 
fitted to high accuracy including small-scale variations. 
For later use in models of actual observations (see sub- 
section \5.2\ , the radial and poloidal velocity field was 
simplified to two or three linear sections depending on 
the particular model. In Figure ^ in grey-scale we show, 
for three different viewing angles, the P-V diagrams of 
the simulation that was also used in Figure (|5|. Note the 
difference in the shape of the first P-V diagram from 
that of the image renderings in Figure (|5| same viewing 
angle). This is a clear indication for deviations from a 
homologous expansion. For a homologous expansion, a 
change in viewing angle should show only a rotating 
structure even in the P-V diagrams. The changing struc- 
ture at different viewing angles confirms the presence of 
deviations. In a reconstruction that incorrectly assumes 
a homologous expansion, the derived 3D structure will 
depend on the viewing angle with deformations along 
the line of sight. 

The contour lines in Figure ^ are those of the simpli- 
fied velocity field model. Together with the reconstructed 
structure in Figure (|5| they show that piecewise linear 
velocity fields can yield an accurate representation of 
the global velocity structure in commonly seen planetary 
nebulae. In the following section we show how this 
result has been applied to a planetary nebula. 

5.2 Reconstruction of the Saturn Nebula 

We exemplify the reconstruction of complex planetary 
nebulae with the case of NGC 7009 |36|, also called 
the Saturn nebula. The top panel of Figure ^ shows 
a color image of the nebula as constructed from three 
narrow-band filter images obtained with the Hubble 
Space Telescope (HST) |1|. In the bottom panel a Shape 
model of the nebula is shown. This is a version of the 
model that includes more structural details than the one 
published in by Steffen et al. |36|. Compared to the 
model in that publication, for the rendering in Figure ^ 
we have added more small-scale features and models of 
additional narrow-band lines. 

To model the object we took into account ground- 
based spectral data of four different spectral lines p9) 
in addition to color HST images. The spectral data are 
composed of 12 position-velocity diagrams rotated at 
intervals of 30° around the central star with a spatial 
resolution of approx. 1 arcsec and a velocity resolution 
of about 7 km/s. 

Sabbadin et al. f29\ used these data to derive the 
structure of NGC 7009 with their tomographic method 
which assumes homologous expansion (see section |zT]). 
Since the spatial coverage of the spectral data is incom- 
plete, they apply an angular interpolation in regions 
of missing data between slits. This produces a unique 
solution for the 3D structure, but at the spatial resolution 
corresponding to ground-based observations, which is 
more than a factor of 5 worse than that of the HST. 



Although the object is not exactly axisymmetric, the 
reconstructed 3D structure of the object shows deviations 
from axisymmetry that are similar to those introduced by 
deviations from a homologous expansion. 

Furthermore, the presence of hot X-ray emitting gas 
within the inner shell |9| leads to the expectation that 
this shell is likely to present deviations from a non- 
homologous expansion. The deviations have been esti- 
mated from a comparison between P-V diagrams and 
images. Although no unique solution was found, the 
distortions of the structure introduced by assuming a 
homologous expansion have been reduced considerably 
by the proposed velocity fields which include a poloidal 
velocity component |32|. They are within the expecta- 
tions from the earlier hydrodynamical simulations. 

In contrast to the direct tomographic reconstruction, 
we interactively modeled the structure as a set of nested 
meshes (Figure ^ and compared the rendered images 
and P-V diagrams with the observed counterparts. The 
interactive flexibility in the modeling of the complex 
structure and velocity field allows it to readily take into 
account complex structures without the need of any 
new coding. Volume meshes have been used for the 
large-scale structures, whereas the thin main shell was 
modeled as a surface mesh. Knots and filaments were 
reproduced using manually added particles on a copy 
of the main shell. The velocity and emissivity distribu- 
tions have been adjusted separately for each mesh using 
piecewise linear functions. 

Observationally, tests of the model can be achieved by 
measuring the expansion component in the plane of the 
sky (which can not be obtained by spectroscopy). Unfor- 
tunately, this requires the detection of the expansion in 
images that have been obtained with a time difference 
of a decade or more with the Hubble Space Telescope. For 
NGC 7009 this has not yet been adequately achieved. 
The velocity vector visualization in Shape provides an 
immediate built-in prediction of the expansion pattern 
in the plane of the sky. Figure (|9| shows the projected 
velocity vectors for one of the models of NGC 7009. 
They clearly show the deviations from a homologous 
expansion, in that the vectors of the main shell and the 
bright symmetric knots do not converge at the position 
of the central star. Future observations of this pattern 
will provide hard evidence in favor or against the model 
published in Steffen et al. | |36) . 

5.3 Nova RS Ophiuchi 

Ribeiro et al. |25 | have used Shape to model the structure 
and expansion of nova RS Ophiuchi after its outburst 
in 2006. In this system a white dwarf star orbits inside 
the outer layers or stellar wind of a red giant star. 
The white dwarf accumulates material from the giant, 
which after some time produces a nuclear explosion on 
its surface. The result is a fast expanding shell around 
the binary system (Figure [lO| top). Nova explosions 
may have speeds of several thousand kilometers per 
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Fig. 7: A research example: The planetary nebula 
NGC 7009 as observed by the Hubble Space Telescope (top) 
and as modeled with Shape by an astronomer (bottom). 
Different colors represent different spectral line images. 
They have been modeled with similar structural meshes, 
but quite different emissivity distributions. This model 
is based on the one published in [[36). 





Fig. 9: The velocity vector components in the plane of 
the sky are shown for the model of NGC 7009 (36). 
The color encodes the velocity component along the 
line of sight, with blue approaching and red receding 
from the observer. While the expansion along the line 
of sight is observable directly using the Doppler-shift 
on spectroscopic emission lines, observing the velocity 
component in the plane of the sky requires the detection 
of position shifts in images that have been obtained with 
a considerable time difference (usually more than five 
years with the HST). The reconstruction of the velocity 
field as shown is a direct prediction for the expected 
expansion pattern of NGC 7009. After its observation it 
can provide evidence for or against the model. 




Fig. 8: In Shape, models are constructed from meshes 
which are generated from simple primitives and de- 
formed using various modifiers. This model of the Saturn 
nebula (NGC 7009) has a number of nested transparent 
shells |36|. 



second, which produces considerable Doppler-shifts in 
the observed spectral lines (Figure [lOj bottom). Some 
emission might then be outside the range of observations 
with narrow-band filters on the Hubble Space Telescope 
and may therefore go undetected. Ribeiro et al. |25| have 
used the spectral rendering filter in Shape to explain such 
''missing'' regions in their HST imaging observations 
of Nova RS Ophiuchi. This feature was introduced in 
the software upon their request. When they do not use 
the HST filter transmission, the object is symmetrically 
double-lobed (Figure 10 top, left). When the filter trans- 
mission is included, one of the lobes largely disappears 
(Figure [lOj top, right). The model with the filter matches 
the observed image much better, although there is still 
some discrepancy in the detailed structure of the larger 



lobe (Figure 10 top, center). The inner bright region 
dominates the ID line profiles (Figure 10 bottom) and 
was very useful to set limits on the orientation of the 
nebula. 

Figure ( pT| shows the mesh of the bipolar nebula that 
they constructed. The inner shaded region was found 
not to expand significantly during the explosion, as 
shown by second epoch observations. The top panels of 
Figure ( pT) show the rendered images on the right and 
left that do and do not include the spectral rendering 
filter, respectively. They can be compared with the ob- 
served image in the middle. The bottom panels compare 
the synthetic spectral line profiles of the final model 
(left) with the observed line profile (noisy line). On the 
right, the range of model line profiles is shown that 
is still compatible with the observations. The difference 
between them is only due to the viewing angle in the 
range from 29° to 40° deviation of the object axis from 
the plane of the sky. 

Using the time modifier with the assumption of ballis- 
tic expansion, Ribeiro et al. |25| have been able to show 
that there is a considerable difference in the expansion of 
the bright inner region and the dimmer bipolar lobes. In 
this developer independent work, the multi-functional 
interactive modeling approach has proved to be espe- 
cially fruitful. 

5.4 Content production for digital media 

Since their emergence, commercial animation systems 
have been applied to produce animated visualizations 
of high spatial and temporal resolution of astronomical 
phenomena in scientific documentaries and feature films. 
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More recently, astrophysical research simulations have 
contributed impressive visualizations of phenomena that 
are impossible to do in commercial animation packages. 
A stunning example is the evolution of a star forming 
region simulated by Henney et al. that was presented as 
part of the show "Journey to the Stars'' by the Hayden 
Planetarium in New York and other digital planetaria. 

Shape has the potential to simulate and visualize a 
variety of astrophysical phenomena for graphical media 
applications, both in animation and stills. A few exam- 
ples of such visualizations can be seen on the Showcase 
page of the Shape website (http://www.astrosen.unam? 
[mx/ shape). With the future development of our software 
and that of computing resources. Shape will become a 
very useful tool to visualize and illustrate astrophysical 
phenomena for a variety of media, including print, tele- 
vision and digital planetaria. 

6 Future Development 

The current and near future development of Shape is 
steered by two driving forces: The first is the type of 
scientific applications that the developers and current 
users are working on. Second, development constantly 
addresses the existing limitations of the software. Lim- 
itations include the availability of processing memory 
in the Java Virtual Machine which mainly translates into 
limited spatial and spectral resolution. Another signif- 
icant limitation is that texture mapping is currently 
not available, which would boost the possibilities to 
model noise structures like complex filaments that are 
observed in many nebular objects. Since the develop- 
ment of the Java version of Shape began, new features 
and changes have been introduced continuously. Current 
developments include explicit radiation transfer for dust 
scattering and spectral lines. 

Long-term plans are guided by potential applications 
that the system has. Such plans include the incorporation 
of interactive hydrodynamical simulations making use of 
multi-processor graphics processing units. Gravitational 
interaction for particle systems in Shape will allow it 
to simulate a variety of phenomena, like interacting 
galaxies and multiple stars. A substantial increase in 
spatial resolution will make the system applicable to 
very realistic modeling and rendering of astrophysical 
objects for educational purposes in planetaria and other 
electronic media. 

7 Conclusion 

We have presented a novel 3D application for the mod- 
eling and reconstruction of astrophysical objects that 
incorporates interactive modeling tools. It considerably 
extends the capabilities of conventional reconstruction 
systems. This is achieved through the support of a 
system of construction and "modifier'' tools that allow 
extremely complex structures and velocity fields to be 
assembled without the need for user programming. A 
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Fig. 10: The top panels show the observed image (center) 
and two different model renderings of nova RS Ophi- 
uchi |25|. The model on the right includes a special 
rendering filter that corresponds to the transmission 
filter on the Hubble Space Telescope. It excludes some of the 
emission with very high velocities along the line of sight. 
The bottom panels show the spectral line profiles from 
the observations (noisy red line) and the final model (left, 
continuous line). The right panel compares the observed 
profile with models at different viewing angles that are 
just consistent with the observations (dashed and dotted 
lines). [Figure reproduced with permission of the authors 
(Ribeiro et al. ||^).] 

number of visualization styles that are common for astro- 
nomical observations can be used to compare the model 
with the observed data. The workflow is enhanced by the 
ability to continuously compare the model to observa- 
tional data during the modeling process as well as by an 
automatized optimization algorithm. The tool has been 
shown to cover the entire modeling and visualization 
pipeline of a common morpho-kinematical modeling 
task, supporting the scientifically accurate reconstruction 
of a wide class of astrophysical objects while keeping a 
convenient and user-friendly interface. The software we 
have presented has been thoroughly tested and applied 
in a number of astronomical research projects, some 
of which we have quoted as examples. In addition to 
scientific research, it may prospectively be applied for 
physically plausible artistic works or for the generation 
of astronomical animations for educational purposes, e.g. 
in digital planetariums. Non-astrophysical uses can also 
be imagined wherever velocity information is observed, 
e.g. in the field of Doppler radar observations of torna- 
dos and other weather phenomena. 

Shape is freely available as a Java WebStart applica- 
tion from its website at http://www.astrosen.unam.mx/] 
shape/ 
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Fig. 11: The outburst of nova RS Ophiuchi has been 
modeled by Ribeiro et al. |25| using an expanding bipo- 
lar nebula (mesh structure) with a nearly inert waist 
(shaded structure). Using the time modifier feature of 
Shape they were able to model the expansion between the 
first and second epoch observations. [Figure reproduced 
with permission of the authors (Ribeiro et al. |25|).] 
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